Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 29 February 2012 (MN MTbX style file v2.2) 



Spectroscopic bulge-disc decomposition: 

a new method to study the evolution of lenticular galaxies 

E.J. Johnston/* A. Aragon- Salamanca/ M.R. Merrifield^ and A.G. Bedregal^'^ 

^School of Physics and Astronomy, University of Nottingham, University Park, Nottingham, NG7 2RD, UK 
^ Departamento de Astrofisica, Facultad de Ciencias Fisicas, Universidad Compiitense de Madrid, 28040 Madrid, Spain 
Institute for Astrophysics, University of Minnesota. 116 Church Street SE, Minneapolis, MN 55455 USA 



29 February 2012 

ABSTRACT 

A new method for spectroscopic bulge-disc decomposition is presented, in which the spatial 
light profile in a two-dimensional spectrum is decomposed wavelength-by-wavelength into 
■ bulge and disc components, allowing separate one-dimensional spectra for each component 

to be constructed. This method has been applied to observations of a sample of nine lenticular 
galaxies in the Fornax Cluster in order to obtain clean high-quality spectra of their individual 
bulge and disc components. So far this decomposition has only been fully successful when 
applied to galaxies with clean light profiles, without contamination from dust lanes etc. This 
has consequently limited the number of galaxies that could be separated into bulge and disc 
components. Lick index stellar population analysis of the component spectra reveals that in 
those galaxies where the bulge and disc could be distinguished, the bulges have systematically 
higher metallicities and younger stellar populations than the discs. This correlation is consis- 
tent with a picture in which SO formation comprises the shutting down of star formation in 
the disc accompanied by a final burst of star formation in the bulge. Similarly, a trend was 
found to exist whereby galaxies with younger stellar populations have higher metallicities. 
The variation in spatial-fit parameters with wavelength also allows us to measure approxi- 
mate colour gradients in the individual components. Such gradients were detected separately 
in both bulges and discs, in the sense that redder light is systematically more centrally con- 
centrated in all components. However, a search for radial variations in the absorption line 
strengths determined for the individual components revealed that, although they can be sensi- 
tively detected where present, they are absent from the vast majority of SO discs and bulges. 
The absence of gradients in line indices for most galaxies implies that the colour gradient 
cannot be attributed to age or metallicity variations, and is therefore most likely associated 
with varying degrees of obscuration by dust. 

Key words: galaxies: elliptical and lenticular - galaxies; evolution - galaxies; formation - 
galaxies; clusters - galaxies; stellar content 



1 INTRODUCTION 

The evolution of galaxies in different environments is still poorly 
understood, as is the relative importance of the different processes 
involved. Lenticular galaxies (SOs) lie between spirals and ellipti- 
cals on the Hubble Sequence, and share many properties with both 
of these morphologies such as the redder colours and older stellar 
populations of ellipticals along with the stellar discs of spirals. As 
a result, SOs are often seen as a transitional phase between spirals 
and ellipticals, and thus understanding their origins is thought to 
be a key stage in understanding galaxy evolution. The many influ- 
ences controlling their evolution are expected to affect the bulges 
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and discs of these galaxies in different ways, and so the individual 
study of these components should provide valuable clues as to their 
evolutionary histories. 

Evidence for the evolution of galaxies along the Hubble Se- 
quence with redshift lies in the morphology-density relation for 
clusters of galaxies. This relation shows that spirals tend to dom- 
inate the less-dense outer r egions, while t he early types tend to 
lie at the cores of clusters l lDresslei|[T980l) . and that the relative 
fraction of spiral galaxies in clusters decreases toward lower red- 
shifts while the fraction of SOs increases ( jPressle r et al. 1997i). 
This finding has been supported by studies using larger samples of 



ataxi es, s uch as Whitmore. Gilmore & Jonel(ll993h. Fasano et alj 
20001) and lDesai etZI l l2007h . Similarly. PWrela et al J | |2004|) have 
also found that late-type spirals are more frequent among iso- 
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lated galaxies, while SOs are more commonly found in denser 
environments such as galaxy groups and clusters. Many scenar- 
ios have been proposed for the transformation of spirals into SOs, 
mostly focusing on the truncation of star formation within the 
galaxy followed by passive evolution into an SO. Processes in- 
clude the removal o f the cold disc gas by ram pressure stripping 
( lGunn&Got^ll972h . the remo val of the hot halo gas by a pro- 
cess usually called starvation l lLarson. Tinslev & Caldwelllll980l : 
iBekki. Shiova & Couchl 2002h. tidal stripping by g a laxy harass- 



ment i Moore et al 199^ 



Moore. Lake & Katj [19981 : iMoore et alj 



Il999h. and star bursts triggered by unequal mass galaxy mergers 
I Mihos & Hemq uist 199 4) and galaxy-clu ster interactions jMerrittl 
ll984lMilleJT986l ; lBvrd & Valtonenlll990l) . 

Since these processes would affect galactic bulges and discs 
differently, many studies have tried looking for clues to SO star 
formation histories by photometric bulge-disc decomposition. The 
paramete rs for the bulge a nd disc can be plotted as the Kormendy 
Relation l lKormendvl[l977l) . where the effective surface brightness 
for each component is plotted against its scale length. This gen- 
erally shows that the surface brightness fall s as the scale length 
increases over a large sample of spirals, but lBarwav et alj | |2009|) 
found that while SOs in the field show this trend, faint cluster SOs 
do not, instead forming a downward scatter where the effective sur- 
face brightness drops significantly at larger scale lengths. This drop 
is thought to be evidence of transformations induced by environ- 
mental effe cts such as minor m ergers, ram pressure stripping or 
harassment. iBarwav et alj ( l2007h also found that in faint SOs the 
bulge effective radius increases with the disc scale length, while 
in brighter SOs it decreases with increasing disc scale length. Sec- 
ular evolution of spirals into lenticulars predicts the former trend 
since the passive fading of the bulge and disc would conserve their 
scale lengths, while the latter correlation would be expected from 
a more turbulent transition. Therefore this variation suggests that 
lenticulars may have evolved in different ways depending on their 
environments. 

By analysing multi-waveband photometry it is possible to 
look at the colour differences and gradients within the bulge 
and disc. Such data can in turn provide information on the star 
formation histories of these galaxies since higher metallicities 
and older stellar populations strengthen the redder light from the 
galaxy. Negative colour grad i ents have been found in t he bulges 
of S Os jTerndrup et al.lll994l ; IPeletier & Balcellslll996l) and spi- 
rals ( lM611enhofdl2004h , showing that the blue light is distributed 
throughout the bulge while redder light is more concentrated in 
the centre. In addition to this it has been found that the discs of 
SOs and spirals are bluer than the bulg es jBothun&Greg3ll99d ; 
IPeletier & Balcellslll993 : lHudson et al.ll2010l) . suggesting that disc 
galaxies have more recent star formation activity at larger radii 
de J ong '1996") or high er metallicities in th eir nuclear regions 
Beckman et al..,1996j : .Pompei & Natalill 19971) . These trends indi- 



cate the presence of age and metallicity gradients across the galax- 
ies, but fail to provide information on whether it represents a gra- 
dient within the individual components, or whether it arises sim- 
ply from the superposition of varying amounts of bulge and disc 
light, where each component contains stellar populations of differ- 
ent ages and metallicities. It also leaves open the possibility that the 
gradients could result from differing amounts of dust extinction at 
different radii, rather than telling us about spatial variations in the 
stellar population. 

To resolve this ambiguity, we have developed a new method to 
apply bulge-disc decomposition to long-slit spectroscopy of nearby 
galaxies. The technique involves dividing the two-dimensional 



spectrum into two one-dimensional spectra, one representing 
purely the bulge light and the other the disc light. This decom- 
position is achieved by fitting the spatial profile of light at each 
wavelength with a disc-plus-bulge model in order to ascertain the 
fraction of the light that should be ascribed to each component at 
that wavelength. By repeating this analysis at every wavelength, 
the individual spectra of bulge and disc can be constructed. These 
"clean" spectra can then be analysed for their ages and metallici- 
ties, as well as any gradients in these quantities, to determine the 
sequence of events in the formation of SO galaxies. 

This paper is intended to present the new method, and give an 
indication of its potential by applying the technique to a relatively 
small sample of SOs from the Fornax Cluster. Section [2] describes 
the Fornax data sample and observations. Section[3]presents the de- 
velopment of the method, and its initial application to the Fornax 
Cluster data. Section|4]analyses these results in terms of the stellar 
populations in the bulges and discs of these galaxies, and Section|5] 
looks at their colour gradients. Section [6] searches for any corre- 
sponding gradients in line indices and hence ages or metallicities. 
Finally, our conclusions are presented in Section[7] 



2 THE DATA SAMPLE 

The initial sample on which we are testing this new method com- 
prises the nine le nticular galaxies in the Fornax Cluster, as classi- 
fied bv lKuntschne r ( 2000), with inclinations above 60 degrees and 
luminosities in the range —22.3 < Mb < —17.3. The data were 
obtained using the 8.2 m AntuA'LT between 2002 October 2 and 
2003 February 24 with the FORS2 instrument in long-slit spec- 
troscopy mode. The slit was aligned with the major axis of each 
galaxy, and was set to 0.5 arcsec wide and 6.8 arcmin in length. For 
all galaxies except NGC 1316, the centre of the galaxy was placed 
half way along the slit in order to obtain information along the en- 
tire length of the major axis; NGC 1316 is the largest galaxy in the 
sample, and it was necessary to offset it along the slit in order to 
cover sufficient sky for background subtraction. The standard res- 
olution collimator was used in the unbinned readout mode, giving 
a spatial scale of 0.125 arcsec pixer\ and the GRIS1400V+18 
grism gave a dispersion of 0.3 1 8 A pixen ^ over a wavelength range 
of 4560 A ^ 5860 A. By analysing arc lines, the spectral reso- 
lution was found to be ~ 4 pixels FWHM, which corresponds to 
a velocity resolution of ^ 73 kms^^ FWHM or a velocity disper- 
sion of ^ 31kms~ ^. Details of the data reduction are given in 
iBedregaletalJ l l2006h . 



3 SPECTROSCOPIC BULGE-DISC DECOMPOSITION 

Conventional photometric one-dimensional bulge-disc decompo- 
sition of galaxies involves measuring the luminosity along the 
galaxy's major axis. The decomposition is obtained by plotting the 
luminosity against radius, and fitting model bulge and disc compo- 
nents to this profile. The disc is generally modeled as an exponen- 
tial. 



Id{R) = lDoeM~R/Ro), 



(1) 



where Ido is the ce ntral surface bri ghtness of the disc and Rq is 
the disc scale length l lFreemanlll970l) . Similarly, the light profile of 
the bulge can be approximated by a de Vaucouleurs profile, 

Ib(R) = iBe exp |-7.669 [{R/Re^^'' - l] | , (2) 



Spectroscopic bulge-disc decomposition 3 




-20 



20 40 60 80 

Radius (arcseconds) 

Figure 1. The light profile of NGC 1375 at 5195 A in black with the best 
fit model in green, compiising the sum of a disc (blue line) and bulge (red 
line) components. 



where i?e is the bul ge effective radius and Ibe is the bulge effective 
surface brightness dde Vaucouleurs|[l953h . For a two-dimensional 
spectrum we can apply exactly the same method to each individual 
wavelength in the data. For example, Fig.[T]shows the light profile 
along the major axis of NGC 1375 at 5195 A, together with the best 
fit achieved using a de Vaucouleurs bulge and exponential disc. 

The only complication here is accounting for the velocity dis- 
persion and radial velocity, where the velocity dispersion of the 
galaxy decreases at larger radii, while the rotational velocity red- 
or blue-shifts the spectra at larger radii relative to that from the 
centre of the galaxy. Therefore, before fitting the light profile at 
each wavelength, the two-dimensional spectrum must be corrected 
such that each spectral feature has the same velocity dispersion and 
radial velocity, thus ensuring that the light profiles for each wave- 
length bin measures the light from the same point on the rest-frame 
spectrum at all radii. Extra care must be taken with this correction 
when a bar is present within the galaxy to account for its kinemat- 
ics. The velocity dispersion was corrected by convolving the spec- 
trum from each spatial location with the appropriate Gaussian to 
bring it up to the maximum value measured within that galaxy. The 
rotational velocity was then corrected by cross correlation, where 
the shift in the wavelength of the spectral features at each location 
was measured relative to those in the peak spectrum. In order to 
minimise the noise in these shift measurements, a rolling average 
was applied to the shifts for each spectrum until the radius where 
the noise dominated the results, at which point the final reliable 
shift was applied to the remaining rows. This gave a smooth ve- 
locity curve from which the necessary shift could be calculated. 
Since the two-dimensional spectra cover the entire length of the 
major axis in all but one galaxy in the sample, the two halves of 
each galaxy spectrum, i.e. the two semi-major axes, were analysed 
separately. This duplication was useful to ensure the results were 
reproducible for each galaxy, and can also provide information on 
spatial asymmetries within the galaxies. 

Once this velocity alignment has been carried out, we can per- 



form the bulge-disc decomposition at all wavelengths. By way of 
illustration. Fig. [2] shows a small section of the models produced 
for NGC 1381 of the bulge, disc and composite galaxy, as derived 
by fitting the spatial profile of the two-dimensional spectrum (also 
shown) wavelength by wavelength. In each profile, the central few 
arcseconds were masked out prior to fitting in order to eliminate 
the effects of seeing, so the peaks of the bulge and disc spectra at 
very small radii in Fig. [2] were not used in the fit. Aditionally, the 
light profiles of each galaxy were also checked by eye to ensure 
that the exponential parts most likely represented the disc rather 
than the spheroid. This was found to be the case in all the galax- 
ies that were decomposed successfully. The results for the bulge 
effective radius and disc scale length for each galaxy from the de- 
composition are given in AppendixiA] Thes e were compared to the 
photometric results of teedregal et alj200^ derived from the Two- 
Micron All-Sky Survey (2MASS) K-band images of these galaxies, 
and were found to be reasonably consistant considering the limited 
model used in this study. 

The total luminosity of the bulge and disc at each wavelength 
can then be calculated by simple integration, using the standard 
results that 



Lfl(A) = 7.227r/e(A)i?^(A) 



for the bulge, and 



(3) 



(4) 



for the disc, where we have now made the dependence on wave- 
length, A, of the various fitted parameters explicit. These quanti- 
ties are simply the modeled spectra of the integrated light from the 
separate bulge and disc components. For example, Fig.[3]shows the 
spectra derived in this way for NGC 1375. The high signal-to-noise 
ratio of the spectra obtained using this integrated light approach is 
immediately apparent, such that one can see that both the H/3 line 
and the magnesium triplet appear stronger in the bulge than in the 
disc; since these features are used as age and metallicity indicators 
respectively, this difference already hints that the bulge contains 
younger stars with a higher metallicity than the disc in this galaxy, 
and we will make a more quantitative assessment of this impression 
in Section|4] 

Application of this technique to the full sample of galaxies 
described in Section|2]revealed the kind of data required to imple- 
ment this analysis successfully. In fact, of the nine galaxies, only 
two (NGC 1381 and NGC 1375) could be reliably decomposed 
into bulge and disc spectra in this manner. A further three galax- 
ies (IC 1963, ESQ 35 8-G006 and ESQ 359-G002) were found by 
iBedregal et"ai] ( |2006|) to have very compact bulges and therefore 
to be disc dominated from very small radii, making it impossible to 
determine a reliable bulge model. For these systems a disc spectrum 
was extracted by assuming that the bulge light was negligible out- 
side of the central masked region, and just fitting a disc component. 
In two further cases (NGC 1380 and NGC 1316) we found that al- 
though the spectra could be decomposed into bulge and disc com- 
ponents, the resulting model did not reproduce the original two- 
dimensional spectrum at all well. Images of these galaxies clearly 
show the presence of major-axis dust lanes in both cases, which 
would affect the light profile, and it has also been suggested by 
ICaon. Capaccioli & D'Onofrid i ll 9941) that NGC 1316 is a merger 
remnant, which would complicate its light profile and therefore 
make it unsuitable for fitting a de Vaucouleurs bulge and expo- 
nential disc. The final two spectra (NGC 1380A and ESQ 358- 
G059) could not be decomposed as their light profiles were again 
too complex for the current simple model. We have therefore as- 
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Figure 2. A section of the two-dimensional spectra for the bulge, disc and composite models, and the original spectrum of NGC 1381, showing how the light 
profile varies with wavelength. The data have been binned such that each step in the radial direction is 12.5", and in the dispersion direction is 6.6A. The 
wavelength range covered is between 5085A and 5243A, increasing towards the right of the diagrams. Note that the peak values of the bulge and composite 
spectra have been reduced by a factor of four in order to allow them to be plotted on the same scale as the original spectrum. 



certained that this approach works best in systems where both 
components are well resolved, and where there is little indication 
of complicating issues like strong dust lanes or extra components 
such as bars. Fortunately, comparison between the original two- 
dimensional spectrum and the model offers a useful a posteriori 
check on the impact of such additional factors. 



4 STELLAR POPULATIONS 

In the integrated light spectrum of a galaxy, the strength of the vari- 
ous absorption lines provides information on the underlying stellar 
population, with hydrogen lines primarily associated with its age, 
while magnesium and iron lines constrain its metallicity. In order 
to obtain quantitative estimates of age and metallicity, these line 
indices can be compared to those predicted by simple stellar popu- 
lation (SSP) models that have been created using stellar libraries of 
the same spectral res olution as the data. In this study, the SSP mod- 
els used are th ose of | Vazdekis et al.l (1201 Ol). w hich uses the MILES 
stellar library dSanchez-Blazguez et alj 20061) . The library spectra 
have a resolution of ~ 58.4 km s~^, and are matched to the spec- 
tral resolution of the data by convolving the stellar spectra with a 
Gaussian of the appropriate dispersion. These data are available as 
a web-based toojjthat allows the user to modify the library spectra 



to the resolution and redshift of their observations, and therefore 
obtain tuned SSP models for their data. 

The line strength i ndices in the d ata were measured using 
the INDEXF software of lCardiell ( l2007h . which uses the Lick/IDS 
index definitions to calculate a pseudo-continuum over each ab- 
sorption fe ature based on the level of the spectrum in ban ds on 
either side jWorthev et allligg^ Iworthev & Ottavianilll997l) . The 
strength of the spectral feature was then measured relative to this 
pseudo-continuum, and the uncertainty estimated from the propa- 
gation of random errors and the effect of uncertainties in the radial 
velocity. 

As one further slight complication, iBedregal et al.l ( 1200 8t) 
found that the galaxies in this sample show traces of emission in 
the H/? line which would reduce the absorption-line strength of this 
line an d skew the results to older ages. To correct for this contami- 
nation, iGonzalezI 119931) identified that the ratio between the equiv- 
alent widths of the [OiiiIasoot and H/3 emissio n features is around 
0.7 in the brightest ellipticals. A later study bv lTrager et al.l ( 120001) 
found the ratio to vary between 0.33 to 1.25 with a mean value of 
0.6. This mean value was used here to estimate the level of emission 
correction necessary for the H/3 index, using 



A(H/?) = 0.6 X EW[Oiii]. 



(5) 



http://miles.iac.es/ 



where the [Oiii]a5007 index was measured from the residual 
spectrum obtained by subtracting the best combinations of stel- 
lar templates from the original galaxy spectrum. These fits were 



Spectroscopic bulge-disc decomposition 5 



NGC1375 



1.4 
1.2 
1.0 
0.8 



3 

E 1-2 

X) 

.2 1.0 



0.8 

0.6 

4500 



Bulge Spectrum 



_ Disc Spectrum 



HP 





5000 



5000 5500 
Wavelength (A) 

Figure 3. The decomposed one-dimensional bulge and disc spectra for 
NGC 1375. 
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Figure 4. Example of the bulge and disc data for NGC 1375 with the 
IVazdekis et al.l ilOlOl) SSP model for the bulge kinematics over-plotted. H/3 
is the age indicator while [MgFe]' is the metallicity indicator The circle 
represents the bulge while the rectangle corresponds to the disc value. The 
error bars represent the statistical uncertainties. 



achieved using the Penalized Pixel Fitting method (PPXF) of 
ICappellari & Emselleii] j2004h . which uses the MILES stellar li- 
brary to model the line-of-sight velocity distribution as a Gaussian 
with a series of Gauss-Hermite polynomials. For this sample, the 
H/3 corrections were never more than 24% of the H/3 index, with 
the majority being less than 10%, so any residuals from the ap- 
proximate nature of this correction are unlikely to compromise the 
results significantly. The PPXF code also used the stellar library to 
measure the line of sight velocity and velocity dispersion for each 
decomposed spectrum, which could then be used to tune the SSP 
models. These values are also given in Appendix IaI 

Having measured and corrected the observed line strengths, 
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Figure 5. Age versus metallicity for the decomposed spectra. Each galaxy is 
represented by a different colour, with circles and rectangles corresponding 
to bulges and discs respectively. Dotted lines join the bulge and disc of the 
same galaxy. The siz e of each symbol reflects the luminosity of the galaxy 
iMadore et ai]ll999h . The filled and open symbols are used to show the 
results from the two semi-major axes in each galaxy, where both could be 
measured. 



SSP models were created for each decomposed spectrum and plot- 
ted using using H/3 and [MgFe]' as the age and metallicity indi- 
cators respectively, where the latter was chosen due to its negli- 
gible dependence on the Q-eleme nt abundance ( lGonzalez|[l993l : 
iThomas. Maraston & Bende3l2003l) . As an example. Fig. |4] shows 
the indices derived for the bulge and disc spectra of NGC 1375, to- 
gether with the grid of predictions for SSP models of differing ages 
and metallicities corrected to the velocity dispersion of this galaxy, 
from which one can read off estimates of the relative age and metal- 
licity for the data. Since these measurements were obtained from 
spectra representing the integrated light over the full bulge and 
disc, they correspond to the global, luminosity-weighted values for 
age and metallicity for each component. The errors shown reflect 
the statistical uncertainties outlined above. Clearly in this case the 
bulge is inferred to be young and metal rich, while the disc is old 
and metal poor. 

Figure|5]shows the results of this analysis applied to all galax- 
ies in the sample. Where both were observed, we show the two 
sides of each galaxy separately, principally as a test of the repro- 
ducibility of the results. Dotted lines join the bulges and discs de- 
rived from the same side of the same galaxy, while points corre- 
sponding to different halves of the major axis of the same galaxy 
share the same colour. In the few cases where the index measure- 
ments lay outside of the SSP model grids, such as in Fig. |4l the 
metallicities were estimated by using an extrapolation method. Due 
to the uncertainties already mentioned, and the fact that different 
SSP models would give different grids, it is important to consider 
the results in Fig. |5]as constraining the relative ages and metallic- 
ities within the data set as opposed to their absolute values. The 
errors shown on the data points correspond to the statistical errors 
given in Fig.|4]plus interpolation errors, while a more realistic mea- 
sure of the uncertainties can be obtained by comparing the results 
for each half of the major axis of each galaxy, as represented by 
closed and filled symbols. It can be seen that both within individual 
galaxies and viewing the data set as a whole, there is clearly a trend 
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Figure 6. Model bulge effective radius and disc scale length plotted against 
wavelength for NGC1375. A strong negative gradient can be seen in the 
disc scale-length, con'esponding to a negative colour gradient across the 
disc, while the gradient in the bulge is much weaker. 



in the sense that the bulges are systematically younger and more 
metal rich than the discs. 



5 COLOUR GRADIENTS 

As well as estimating the global ages and metallicities of the bulges 
and discs, this spectral decomposition method also provides an ap- 
proximate measure of any gradients in the properties of compo- 
nents. Specifically, if there were a colour gradient in a component, 
such that, say, red light was more centrally concentrated than blue 
light, then the characteristic size-scale of the component would be 
smaller in red light than in blue light. Because we determine the 
length-scales of both disc and bulge as a function of wavelength, 
it is straightforward to see whether any such gradients exist. For 
example. Fig. |6] shows the bulge and disc length-scales across the 
range of wavelengths observed for NGC 1375. In both cases, there 
is a negative gradient, indicating that the centre of each component 
is redder than its outer parts, although the effect is clearly stronger 
in the disc than the bulge. 

In order to quantify this effect in terms of more conventional 
colour gradients, we can calculate the bulge effective radius and 
the disc scale-length at the the central wavelengths of the B and V 
band filters f rom the Johns on-Cousins system [4450 A and 5510 A 
respectively ( lBesselliri990h l. For the bulge, the ratio of light in the 
B-band to that in the T/-band is given by 
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Figure 7. The colour gradients in the bulges (circles) and discs (rectan- 
gles) of each galaxy against the galaxy's absolute S-band magnitude. The 
colours of the points are as in Fig.|5] 



Differentiating with respect to i?^'''', we find 



d{B - V) 



= 8.3 



R. 



1/4 



1/4 



(8) 



This equation can then be integrated out from the centre to any 
radius, to obtain a change in colour. 



(B - V)r - (B - V)o = 8.3 



pl/4 



R 



1/4 



(9) 



As a characteristic radius, we set R = R^v in Equation and 
then divide by log Rcv to create an appropriate gradient quantity. 

We can also obtain a set of analogous equations in order to 
define a gradient in the disc: 



Ib_ ^ (Iob\ 
Iv \Iov J 



exp 



1 



R 



1 

Rov 

loB 



(10) 



B -V = const - 2.5 logi 



d{B - V) 
d{R) 



= 1.09 



1 

Rob 



Rov ) ' 



(11) 



(12) 
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const — 2.51ogio ( 7^ 
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1/4 



R, 



1/4 



1 



1/4-1 



R 



1/4 



(6) 



(7) 



(B - V)r - (B - V)o = 1.09 



Rob Rov 



R. 



(13) 



The colour gradients obtained from these formulae using the 
observed variation in scale-length with wavelength for each com- 
ponent of each galaxy are shown in Fig. |7] In all cases, the 
colour gradients are negative, indicating that the centres are red- 
der than the outskirts, in both bulge and disc components. There 
also appears to be a trend that the gradients are stronger in fainter 
galaxies. The range of values obtained are directly comparable 
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Figure 8. The line gradients for the bulge and disc in each galaxy, calculated using AIndex*/A logj^Q(R). 



to th ose found in previous studies of early-type galaxies, with 
iLa Ba rbera & de Carvalho ( 2 0091) reporting a typical gradi ent of 
~ -0.071 ± 0.003. iRoche. Bemardi & Hvdd l llOlOh find- 
ing a range of -0.8 < < 0.0, and lSuh et al.l feOlOl) report- 
ing values of -0.4 < ^gf^j < 0.0. Iden Brok et Zl l l201lh simi- 
larly detected gradients within the range —0.2 < ^^^"^ < 0.0 for 
a sample of early-type galaxies, but found that SOs show weaker 
gradients than ellipticals, suggesting that the underlying mecha- 
nisms may differ. 

Attempts have also been made to ascertain the c a uses o f such 
gradients. For example, iLa Barbera & de Carvalhd ( |2009|) con- 
cluded that the main contributors to the colour gradients were the 
radial variations in metallicity over the galaxies, and that, while a 
small positive age gradient was also present, its contribution to the 
colour gradient was negligible in comparison. The spectral decom- 
position method presented here provides us with a new mechanism 
for seeing directly in the individual components whether these gra- 
dients are associated with changes in age or metallicity, and we now 
describe how this information can be exploited. 



6 LINE INDEX GRADIENTS 

If there were a gradient in the strength of a particular absorption line 
within a single component, one would expect that the scale-length 



as determined within the absorption feature would differ from that 
of the surrounding continuum, which, in turn, would show up in 
the spectrum of that component generated from the best-fit model, 
like those shown in Fig.|2] We have therefore analysed these model 
spectra by calculating the Lick indices from them as a function of 
radius. A logarithmic gradient was then generated for each index 
by calculating the magnitude version of the indices. 

Index- = -2.5 log,o(^l-^^j, (14) 

and measuring its variation with A log]^Q(i?). The resulting gradi- 
ents for the bulge and disc of each galaxy are presented in Fig.[8] 

Since any signal in this plot arises from a variation in the 
component's scale-length at the wavelength of the index, the error 
on each measurement was assessed by calculating the indices that 
one obtains by treating the plots of scale-length versus wavelength 
(such as those shown in Fig.[6ll as spectra, and extracting Lick in- 
dices from them: if there were a signal, then the value of scale- 
length in the central band should differ significantly from the value 
in the pseudo-continuum, leading to a non-zero index. To measure 
the significance of any such measurement, and hence the appro- 
priate size of the error bars in Fig. [8] the indices measured at the 
true index wavelengths in plots like Fig. |6] were compared to those 
obtained from random locations in these plots. These results were 



8 E.J. Johnston et al. 



ESO 359-GG02 





1 


4 
















1 


2 








OJ 








1 











c 






u 





8 


o 




z 







o 

% 10 

O 
1-1 

i 8 
6 



Disc Spectrum 



Mg triplet 



4 t 

4500 



6000 



5000 5500 
Wavelength (A) 

Figure 9. The decomposed disc spectmm for ESO 359-G002 (top) and a 
plot of its disc scale length against wavelength (bottom). A feature is present 
at the wavelength of the H/3 line on the scale-length plot, implying that there 
is a gradient in its line strength. 



found to be consistant with those o f lBedregalet"ai] ( [20T l|), in which 
the line strength gradients were measured over the whole galaxy. 

As can be seen from Fig.IH the H/3 index of ESO 359-G002 
shows a gradient detected with a > 4ct significance, which would 
suggest the presence of an age variation across its disc. The origins 
of this detection can be seen in Fig.|9] where the scale length mea- 
sured at the wavelength of H/3 differs from that in the surrounding 
continuum, just as described above. The effect of this variation on 
the reconstructed spectra is illustrated in Fig.fTO] which shows the 
derived H/3 feature at radii of 2 arcsec (the inner seeing limit of 
the galaxy) and at 8 arcsec (the disc scale length). The SSP model 
for this galaxy was used to translate the H/3 line index gradient 
into an approximate age gradient of A logjQ(age)/A logjQ(i?) ~ 
1.1 ± 0.2. This corresponds to ages of 1.86 ± 0.14 Gyrs and 
8.77 ± 0.16 Gyrs at the respective inner and outer radii shown in 
Fig. [To] Note, however, that such an age gradient is in the wrong 
sense to explain the colour gradient seen in this system. 

In trying to explain the observed colour gradients, the situa- 
tion is no more helpful in all the other galaxies. As is evident from 
Fig. [8] none of the components in any of these systems display 
significant gradients in any of their indices. We can translate these 
limits into limits on the variation in the properties of the stellar 
population between the centre and each component's characteristic 
radius {Rc for bulges and Ro for discs) as above. Typically, this 
analysis yields an upper limit of a ~ 30% change in age and metal- 
licity between these locations. This correspon ds to a colour gra- 
dient of I '^l:'^zV I ~ 0.02 using the models o f iBruzual & Charloll 



J 1 I ^ log IX I 

( l2003h . which cannot explain the sizable colour gradients found in 
Section |5] It would therefore appear that some other factor, most 
likely dust reddening that varies with radius, must be the underly- 
ing cause of the colour gradients. 



7 DISCUSSION 

We have presented a new method for analysing SO galaxy spectra 
by decomposing their major-axis light into bulge and disc compo- 
nents on a wavelength-by-wavelength basis, in order to construct 
clean, high-quality spectra for each individual component. Applica- 
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Figure 10. The H/3 feature in the reconstmcted spectra of the disc of 
ESO 359-G002 at the radii of the inner seeing limit (black), and the disc 
scale length (red). 



tion of this method to a preliminary sample of data from the Fornax 
Cluster has revealed that as long as the galaxies are well-described 
by this two component model, and the bulge region is sufficiently 
resolved to allow a reliable structural fit, it is quite possible to spec- 
trally decompose galaxies in this way. 

Even this initial small sample reveals some clear systematic 
trends in the properties of Fornax SOs. Specifically, analysis of the 
Lick indices of the individual components shows that in the two 
galaxies that could be decomposed into two components, the bulges 
are systematically younger and more metal rich than the discs. 
This is also found to be the case more globally when comparing 
galaxies across the sample. Since these quantities are luminosity- 
weighted averages for the individual components, we cannot rule 
out the probability that the differences arise from relatively small 
amounts of late star formation, but they nonetheless point to real 
diffeences in star formation history. The positive age gradient be- 
tween the bulge region a nd disc region has pre viously been docu- 
mented in this data set bv lBedregal et al 1 201 ll). and has been see n 
in other lenticu la r galax ie s bvlFisher, Franx & IllingworthI (fl99^. 
iBell & de Jo"n3 ||2000|). iKuntschneJ (l2000h, [Mac Arthur et al J 
1 20041) and lProchaska Chamberlain et al.l ( 1201 ih . but the new anal- 
ysis technique allows us for the first time to cleanly measure the 
relative ages of the individual components, and also to optimise the 
signal-to-noise ratio of the data, and hence the precision of the mea- 
surements, by combining all of the two-dimensional spectral data 
into just two one-dimensional spectra. 

By analysing the variation in the bulge and disc components' 
scale-lengths with wavelength, we were also able to make approx- 
imate measurements of the colour gradients in the individual com- 
ponents, which turned out to be systematically negative, with the 
centre of each component redder than its outskirts. Again, such gra- 
dients have been seen before, but here for the first time we can show 
clearly that they exist in the individual components, not just arising 
from the superposition of monochromatic individual components 
whose relative contributions vary with radius. 

By analysing any variations in the best-fit structural parame- 
ters at the wavelengths of the various Lick indices, we have fur- 
ther been able to measure gradients in age and metallicity. In only 
one galaxy was a significant age gradient detected, which interest- 
ingly was found to be in the wrong sense to explain the system's 
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colour gradient. In all other components of all other galaxies, no 
discernable gradients were detected. The upper limits these mea- 
surements imposed were also sufficient to rule out the possibility 
that any of the observed colour gradients could be due to variations 
in the stellar population. It therefore seems likely that the red cen- 
tres of all components must be attributed to centrally-concentrated 
dust in these systems. 

Although the first results presented here already illustrate the 
potential power of spectroscopic bulge-disc decomposition, there 
are still plenty of areas for further development. One could extend 
the model by, for example, introducing a more general Sersic pro- 
file instead of the de Vaucouleurs law to fit a wider range of bulge 
types, or one could introduce to the model additional components 
such as a bar, in order to see where the stellar population proper- 
ties of such components fit into the developing picture of SO galaxy 
formation. 

It would also be helpful to apply this analysis to a somewhat 
larger sample obtained in a different cluster, both to improve the 
statistics and to search for variations between SO properties in dif- 
ferent locations. To this end, we are currently analysing long-slit 
spectra from a further sample of 2 1 SO galaxies in the Virgo Clus- 
ter, and look forward to presenting the results shortly. 

In addition, the method could be straight-forwardly extended 
to IFU data, where the effects of small scale variations, such as 
those due to dust lanes, can be reduced, and the star formation 
history throughout the two-dimensional structure can be measured 
more robustly. We are currently modifying this technique in order 
to apply it to such data. 
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APPENDIX A: TABLES 

In this appendix we include tables with different parameters 
for each galaxy in the sample. Only NGC 1381 and NGC 1375 
could be decomposed into bulge and disc components, IC 1963, 
ESO 358-G006 and ESQ 359-G002 gave disc only fits, and 
NGC 1316, NGC 1380, NGC 1380A and ESO 358-G059 could 
not be fitted with the current model. The errors are presented in 
parentheses next to the values, and represent one sigma uncertain- 
ties in the measurements. The results for J?e, Ro, Vlos and g are 
given for both halves of the major axis in each case, where Vlos 
and a are the line-of-sight velocity and velocity dispersions after 
the corrections described in Section[3] 





Ma 


Re 


no 










[arc sec] 


[arcsec] 


[kms~^] 


[kms^^] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 1381 


-19.1 (0.1) 


6.7 (0.2) 


25.1 (0.2) 


150.5 (12.6) 


1718.5 (8.8) 






7.5 (0.2) 


26.4 (0.3) 


157.9 (14.0) 


1715.0 (9.0) 


NGC 1375 


-18.2(0.2) 


5.0 (0.4) 


16.4 (0.3) 


115.0 (9.9) 


785.2 (6.6) 






3.5 (0.4) 


16.8 (0.5) 


116.5 (9.4) 


778.3 (6.5) 


IC1963 


-18.5 (0.1) 




15.5 (0.3) 


59.0 (2.3) 


1635.5 (1.7) 








15.0 (0.3) 


61.8(1.8) 


1631.7(1.6) 


ESO 358-G006 


-17.5 (0.2) 




10.5 (0.2) 


98.8 (5.7) 


1264.2 (4.3) 








10.4 (0.2) 


99.6 (6.4) 


1258.5 (4.7) 


ESO 359-G002 


-17.3 (0.2) 




7.3 (0.5) 


58.9 (3.3) 


1439.7 (2.5) 








7.4 (0.5) 


62.7 (3.2) 


1457.3 (2.6) 


NGC 1316 


-22.3 (0.1) 










NGC 1380 


-20.6 (0.1) 










NGC 1380A 


-18.1 (0.2) 










ESO 358-G059 


-17.4(0.2) 











Note. Column (2): absolute magnitude in B-band according to lMadore et alj ( Il999l) : Column (3): bulge effective radius, 
measured at the central wavelength of the V-band filter (5510 A); Column (4): disc scale length, measured at the central 
wavelength of the V-band filter (5510 A); Column (5): velocity dispersion; Column (6): line-of-sight velocity. 



